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Abstract 

We present final measurements of tliirteen charmless hadronic B decay modes from the CLEO 
experiment. The decay modes include the ten tttt, Kit, and KK final states and new limits on 
dibaryonic final states, pp, pA, and AA, as well as a new determination of the ratio B{B 
DK)/B{B Dtt). The results are based on the full CLEO II and CLEO III data samples 
totalling 15.3 fb~^ at the T(4S'), and supercede previously published results. 
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I. INTRODUCTION 



Charmless decays of B mesons may proceed by 6 — > 6 — >■ s, or 6 ^ d transitions. 
The latter two mechanisms require flavor changing neutral currents which are not present at 
tree level in the Standard Model, and therefore must occur through higher order processes 
such as the penguin mechanism. Such processes involve loops, which can open the window 
for particles and physics outside the Standard Model. Even in the absence of such new 
physics, interference among competing amplitudes for a given decay mode can be exploited to 
measure CKM phases. There is a significant body of literature 0] on the use of the branching 
ratio B{B — > K~tt~^) and other charmless modes to determine or constrain the CKM angle 
7, the phase of Vub in conventional representations of the CKM matrix. Compared to the 
methods of extracting 7 that are based the B DK decay modes|2], these approaches based 
on charmless decay modes are less clean theoretically, but more promising experimentally 
because the event yields are significantly higher. 

Recent work on two-body charmless decay modes suggests that the unitarity triangle 
may be constructable entirely from charmless modes, without recourse to the traditional 
constraints involving B mixing measurements, CP asymmetry m. B ^ J/ijjKg, or CP viola- 
tion in kaon decays. The charmless modes therefore offer an independent approach to probe 
CP violating effects in heavy quark decay. Significant disagreement between these two ap- 
proaches, if found in experimental results, would directly challenge the Standard Model and 
its fundamental statement that all CP violating phenomena stem from a single phase in the 
CKM matrix. Early results based on current data are already available|3], and do indicate a 
degree of inconsistency. In this paper we present new experimental data on charmless modes 
and note that these data enhance rather than ameliorate the discrepancy. 

CLEO has previously published several papers j^i] reporting measurements of charmless 
hadronic B meson decay modes, including searches for charmless baryonic final states, with 
the data of the CLEO II experiment. Here we report corresponding measurements in the 
new CLEO III data with results for three vrvr modes, B — * vr+vr", tc^tc^, tt^tt'^, four Kir 
modes, B K+n-, K+n^, K°7r+, KOtt^ three KK modes, B K+R-, K^R-, K^W, 
and three dibaryonic modes, B pp, pA, AA. We also merge CLEO II and CLEO III 
results to determine a final measurement for each mode based on the full CLEO data set, 
which hereby supercedes our previous publications. Recent measurements from BABAR 
and Belle are in excellent agreement with ours We also report a new measurement of the 
ratio of branching ratios, B{B — >■ D^K~)/B{B — >■ D°7r~). 

Here and throughout this paper charge conjugate modes are implied. We also make use 
of the notation to represent a charged hadron that may be either a kaon or pion. 



II. THE CLEO DETECTOR AND DATASETS 

CLEO is a general purpose solenoidal magnet detector operating at the Cornell Electron 
Storage Ring (CESR). The latter is a symmetric-energy storage ring tuned for the data 
sets discussed here to provide center of mass energies near the T(4S). At y/s = ^■Y'(^4g^ 

the hadronic cross section is approximately 4 nb, with 1 nb of e^e^ T(4S) BB 
and 3 nb of four-flavor continuum e~^e~ —>■ qq. In the CLEO HI running period, July 2000 
through June 2001, we obtained an integrated luminosity of 6.18 fb~^ at the T(4S) and 2.24 
fb^^ off-resonance, i.e., just below the BB threshold. The off-resonance data are used for 
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background determinations. The on-resonance data corresponds to N^g = (5.73±0.47)xl0^ 
T(4S) decays. The corresponding numbers for the CLEO II running period (1990-1999) are 
9.13 fb^^ ((9.66 ± 0.19)xlO^T(4S) decays) and 4.35 fb'^ Differences in the Nbb yield per 
unit integrated luminosity reflect differences in run conditions. 

The CLEO III detector [i] differs from the CLEO II detector most notably in the 
inclusion of a ring-imaging Cherenkov device (RICH) [8] which provides particle identification 
at all momenta above the Cherenkov threshold. Even at the highest momenta relevant for 
B physics, about 2.8 GeV, the RICH separates kaons and pions by 2.3 standard deviations. 
Measurements of specific ionization {dE/dx) in the drift chamber provide an additional 
2.0 standard deviation separation at the highest momenta. Charged particle tracking is 
done by the 47-layer drift chamber and a four-layer silicon tracker which reside in a 1.5T 
solenoidal magnetic field and provide momentum resolution described by (cXp/p)^ = A"^ + 
B'^p^ with A ^ 0.005 and B ^ 0.001 GeV~^. The absolute momentum calibration is 
confirmed by comparing the invariant mass of standard decays J/ip ^ f^'^fJ'^j ~^ K^n^ 
with PDG values@. Photons are detected using a 7800-crystal CsI(Tl) electromagnetic 
calorimeterwhich is unchanged between CLEO II and CLEO HI. 



III. ELEMENTS OF THE ANALYSIS 

The T(4S) is produced at rest in the lab frame and decays with low Q value to a pair of 
B mesons that travel non-relativistically, with pr ~ 300 MeV. In this analysis we assume 
equal rates of B^B~ and B^B^ production jlCI]. All decay modes studied in this paper 
are two-body or quasi-two-body modes. Apart from the modest ±150 MeV doppler shifts 
resulting from the motion of the B mesons in the lab frame, the daughter particles are 
nearly monochromatic, and, up to resolution smearing, jointly carry the full beam energy 
El, and have invariant mass equal to the B mass Mb- The approximate monochromaticity 
of the daughters simplifies particle identification and energy resolution, and helps keep the 
associated systematic errors low. The other B in the event decays into, on average, five 
charged and five neutral particles, distributed uniformly in the detector acceptance. The 
principal background to the analysis comes from the non-6 hadronic data, e^e~ qq, with 
qq = uu, dd, ss, and cc. High momentum, back-to-back particles are typical in such events, 
and some have invariant masses and total energies close to or in the signal region of the B 
events. Fortunately, distinctive event topologies separate most of these background events 
from the signal. 

This analysis has two principal parts: (a) the application of hard selection criteria to 
obtain signal-like events, based on kinematics, event topology, and particle identification; 
(b) the application of an unbinned extended maximum likelihood fit to the surviving event 
ensembles to extract the yields of signal and background(s) for each mode. The likelihood 
fit allows us to make maximum use of available information, while avoiding efficiency losses 
that further selection criteria would entail. 

For the purposes of reconstruction, the CLEO HI dataset reported on here is divided into 
two subsets of roughly equal integrated luminosity, which we will call Set A and Set B. The 
distinction has ultimately no significant effect on the results, but because of changes in event 
reconstruction algorithms between the two sets, there are slight differences in resolutions 
and efficiencies - mostly affecting modes with charged particles - that we treat separately 
until the final CLEO HI results are reassembled at the end. We provide in Table |1] some 
informative comparisons between Set A and Set B. 
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TABLE I: Features of Set A and Set B. 



OllPlTlfiifiV 


Set A 


Set B 


Fraction of total -/Vg^ 


00 /o 


4o/o 


Track Resolution 






'A' Coefficient 


0.0055 


0.0044 


-D uoemciem; (jG v j 


n nni i 

U.UUii 


n nni n 

U.UUiU 


J\ TV iViOCLc 








Z. 1 


z. / 


(7^^(MeV) 


22 


19 


Efficiency 


oo/o 


40/0 


J\ TT iviocie 






(MeV) 


o.i 


Q 1 
O.i 


aA^(MeV) 


31 


31 


Efficiency 


33% 


35% 


vr^vrO Mode'^ 






^Mb (MeV) 
c7^^(MeV) 


3.6 


3.6 


43 


43 


Efficiency 


22% 


22% 



"Resolutions are given as average of low-side and high-side half-resolutions. 

Charged track and photon candidates are required to satisfy loose quahty requirements 
which reject poorly determined candidates while retaining high efficiency for real tracks and 
showers. Kg candidates are selected from pairs of charged tracks forming well-measured dis- 
placed vertices with a TT~^7r~ invariant mass within three standard deviations of the nominal 
Kg mass. In addition the vertex must satisfy Iry^xl > 5mm in the transverse plane, and 
Pj^o -tytx > 0. The Kg —>■ tt^tt^ mode is not used. A candidates consist of pvr pairs with in- 
variant mass within three standard deviations of the nominal A mass. Pairs of photons with 
an invariant mass within 2.5 standard deviations of the nominal 7r° mass are kinematically 
fit with the mass constrained to the nominal 7r° mass. 

A. General Event Selection 

Candidates for rare B decay events are selected for further analysis on the basis of two 
kinematic variables and one event-shape variable. For each candidate, we construct the beam 

constrained B candidate mass Mb = \J {El — p ^) where is the beam energy, and p is the 
momentum of the candidate computed from the vector sum of the daughter momenta. For 
real B mesons \p\ <C Ei, and the width of this distribution is dominated by the ~ 2.5 MeV 
intrinsic beam energy spread. The beam energy is determined run by run from CESR lattice 
information, and slight corrections are applied afterward to ensure that the observed B~ 
mass in B~ — >■ D^7i~ events matches the accepted valueQ]. In addition we compute the 
energy balance variable AE = E — Ej, where E is the sum of the daughter energies. The 
width of this distribution is about 20 MeV in all charged modes, as determined by the 
momentum resolution of the tracking systems, and is about 40 MeV in modes involving 
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neutral pions. 

Any candidate with \AE\ < 400 MeV and Mb > 5.2 GeV is kept. An additional re- 
quirement on cos^sph, the cosine of the angle between the sphericity axis of the candidate 
and the sphericity axis of the rest of the event 0], is used to reject the dominant e~^e~ qq 
background. All candidates must satisfy the requirement |cos6'sph| < 0.8, which rejects 
approximately 80% of the background while retaining nearly 80% of the signal. 

B. Particle Identification Requirements (CLEO III) 

In the case of a candidate mode involving one or more charged pions or kaons, such as 
B — > Ktc ot B ^ TTTT^, each charged track must be positively identified as i^" or vr. The 
pattern of Cherenkov photon hits in the RICH detector is fit to both a kaon and pion 
hypothesis, each with its own likelihood Ck and Cj^. The mean number of photon hits 
entering the fit is twelve, and we require a minimum of four. Calibrated dE/dx information 
from the drift chamber is used to compute a fo^^ kaon and pion hypotheses. The RICH 
and dE/dx results are combined to form an effective difference, 

AKn = -2lnCK + 2lnC^ + X^K-x\- (1) 

Kaons are identified by Ak-k < and pions by —Axn < ^tt, with values of Sk and 
chosen to yield (90 ±3)% efficiency as determined in an independent study of tagged kaons 
and pions obtained from the decay D*~^ —>■ tt^D^ {D^ — > K~n~^). With this choice of 6k 
and Jjr; the misidentification rate for kaons faking pions (pions faking kaons) is 11%(8%) at 
momenta around 2.6 GeV. 

For candidate modes involving protons, positive proton identification is required. dE/dx 
does not distinguish well between protons and kaons at the ~ 2.5 GeV momenta of interest, 
however, so the proton-kaon separation is achieved with a discriminant based only on RICH 
information: Axp = — 21n£p + 2\n Ck < 6p. In this case 6p is chosen to yield proton 
(antiproton) identification efficiency of 76 ± 1(72 ± 1)% with a kaon fake rate of 1%, as 
determined in an independent study using tagged kaons from the D*~^ sample as above, and 
protons from A — pvr decays. 

C. Event Selection for CLEO II modes 

We present three results for which we also analyzed the full CLEO II data set, namely 
B — s> K^K^, B — s> AA, and B pA. The selection required that the Kg vertex is 
separated from the beam spot by more that 3 sigma (5.5 sigma for CLEO II. V for which the 
innermost drift chamber was replaced with a 3 layer silicon vertex detector). The candidate 
mass must lie within 10 MeV of the nominal Kg mass. We require that the Kg flight 
direction points to within 3 sigma of the beamspot . 

The protons in the pA final state must be compatible within 3 sigma with a proton 
dE/dx hypothesis and incompatible with both the electron (calculated from calorimeter 
information) and muon (calculated from muon chamber information) hypotheses. We require 
that the A candidate mass lie within 10 MeV of the nominal A mass, the vertex be at least 
5 mm removed radially from the beam spot, and the of the vertex fit be less than nine. 
There is no particle identification applied to daughter particles of the A decays. 
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IV. ANALYSIS VARIABLES 



Events which meet all the requirements described in the preceding paragraphs are now 
used in a likelihood fit to extract signal yield. We characterize each candidate event by four 
variables: the mass and energy variables introduced above, Mb and AE, the flight direction 
of the candidate B, and a Fisher discriminant flT] . 

The flight direction is given by cos 6^ = p ■ z where p is the vector sum of the daughter 
momenta and z is the beam axis. Since the vector T(4S) is produced in e~^e~ annihilation 
it has a polarization = ±1, and the subsequent flight direction of the pseudoscalar B 
mesons is distributed as \Y^'^{9, (f))\'^ ~ sin^^ = 1 — cos^ ^. Background events are flat in 
this variable. 

The Fisher discriminant is used to refine the separation of signal and e"^e^ — > qq back- 
ground that is initially addressed by the hard cut on cos ^sph in the general event selection. 
The Fisher discriminant, Xjr, is a linear combination of fourteen variables, with coefficients 
chosen to maximize the separation of signal and background events. The optimization pro- 
cedure uses Monte Carlo events for the signal and off-resonance and {Mb,AE) sideband 
data events for the background. As in previous CLEO publications 0] the component vari- 
ables include the direction of the thrust axis of the candidate with respect to the beam 
axis, cos^thr, and the nine conical bins of a "Virtual Calorimeter" whose axis is aligned 
with the candidate thrust axis. A fuller description of the Virtual Calorimeter is available 
in a previous publication l3| . Note that cos 6'thr and cos 9b are quite different quantities. 
For two body decay B — »• XY, one has simple closed form expressions: cos 6b = p ■ z with 
p = px + Py, whereas cos 6'thr = q • z with q = px — Py- 

In addition, we take advantage of the high quality particle identification in CLEO III 
to augment the Fisher discriminant with information on the presence of electrons, muons, 
protons, and kaons in the event. The momentum of the highest momentum electron, muon, 
kaon, and proton are used as inputs to the Fisher discriminant. For these purposes we 
need only rudimentary particle identification criteria. If any of the possible particle type 
hypotheses has no corresponding track identified (which is very often the case), a value zero 
is used as the input to the Fisher discriminant. 

The Fisher variable thus defined provides discrimination between charmless B decay 
signal modes and qq background at a level equivalent to two gaussian distributions separated 
by 1.4(7, and is independent of the details of the signal mode for all the modes studied here. 



V. LIKELIHOOD FIT 



A. Fit Components 

With the four analysis variables Mb, AE, xjr, and cos^b? we characterize each event 
in terms of normalized probability distribution functions {PDFs): M.^k{,Mb), £^xn{AE), 
^fxnixyr), and C^k{cos9b). The thirteen different charmless decay modes to be fit will 
in general have contributions from (a) signal, (b) qq background, and (c) cross-feed from 
other B modes. Subscripts /x and k identify the particular decay mode (/i) and the type of 
contribution (k). The probability that a given event characterized by {Mb, AE, xjr, cos^b) 
is an event of component type k of decay mode /z is then given by the product of PDFs 

PfxK = M^n{MB)S^n{^E)J^^^{x^)C^t,{coseB). (2) 
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We determine the yields n/^ of signal, qq background, and cross-feed background in decay 
mode /i by maximizing the extended likelihood function with respect to the yields n/^: 

C.^i{nsig,nqq,n^^f^ed) = ]J i^nt^P^A. (3) 

K events \ K / 

The qq background is the dominant background source in all cases, and in only five of the 
fifteen modes do we need to include any cross-feed backgrounds. Four of these are due to the 
~ 10% K/iT misidentification probability. In fitting B — > tt^tt' and B ti^tt'^ we include 
components for B — > K^t^^ and B K^tt^, respectively; in i? — K^K^ we include a 
component for B — > i^+7r~; and in i? — > D^K~ we include a component for B D^n^. 
Although the cross-feed backgrounds arise from mistaken particle identification, they are still 
distinguishable from the signal through S{AE), which is shifted by about 50 MeV relative 
to the signal PDF. The cross-feed fits are only for background removal and the yields are 
not used in any other signal determination. 

The fifth mode requiring a cross-feed component is B^ tt^ti^. In this case a small 
contribution from B^ p^n^ arises when the charged pion has very little momentum in 
the lab frame. Although the missing particle also shifts AE by at least one pion mass, 
resolution smearing leaves a small tail in the signal region. The treatment here is the same 
as in our previous publication on B^ 7r°7r°jl|. We note also that potential feedthrough of 
B^ —>■ p^TT^ into B^ n^TT^ is smaller than in the tt^h^ case because the low-side resolution 
smearing is less for the vr+vr^ mode, and because the ratio B{B^ — > '7r"'"7r°)/i3(-B"'" — >■ p"'"7r°) is 
larger than B{B~^ — >■ n^7r^)/B{B^ — >■ p+7r°). Monte Carlo studies confirm these observations 
and we therefore do not include this term in the 7r"'"7r° fit. 



B. PDFs 

We parametrize the PDFs with various functions and combinations of functions which 
are listed below. In each case the parameters of the function are determined from a fit to 
signal Monte Carlo event samples for the signal component and cross-feed component (if 
there is one), and from a fit to off-resonance data for the qq background component. These 
parameters are then fixed for all subsequent fitting procedures so the only free variables in the 
likelihood function Eq. |21are the signal and background yields. There is of course underlying 
uncertainty in the parameter values which fix the PDF shapes, but this uncertainty is 
systematic in nature and will be discussed later in section IVII All functions are normalized 
to unit area over the accepted range of the free variable. 

• Gaussian (G): used for Mb and AE signal component PDFs that do not involve 
neutral pions. The parameters are the mean and width. 

• Asymmetric Gaussian (Q): used for Mb and AE in modes where neutral pions appear. 
Fluctuations in the measured 7r° energy are intrinsically asymmetric - with a longer 
tail on the low energy side - because of energy leakage out the back of the Csl crytals 
in the electromagnetic calorimeter. The parameters are the mean and separate left 
and right widths. 

• Linear (L): used for qq backgrounds in AE. The free parameter is the slope. 
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• ARGUS (A): used to characterize the Mb shape of qq backgrounds [l^]- A(Mb) = 

y^l — Ms^/-E| exp(— A(l — Mb"^ /El)). The parameter A governs the turn-over of the 
shape and the slope at low values of Mb- The beam energy E'f, determines the endpoint 
of the qq Mb spectrum. Over the course of CLEO III data taking this end point 
clusters around several close but not identical values. In practice we form a sum of 
three ARGUS functions with different endpoint values, weighted by the corresponding 
integrated luminosities. In addition, to account for run-to-run beam energy variation, 
we convolve each ARGUS function with a Gaussian of width a ~ 0.7 MeV in Eb. 

• Breit-Wigner (7^): used in Fisher parametrizations to describe non-Gaussian tails. 
Parameters are mean and width. 

• Fisher {Fq): a linear combination of functions used to characterize the qq background 
Fisher shape. It is primarily an asymmetric Gaussian (87% of the area), but includes 
an additional Breit-Wigner (9%) with the same mean, and a small symmetric Gaussian 
(4%). 

Table im lists the PDFs used for each fit component of each mode. The fourth fit variable, 
cos^^B is in all cases taken to have the functional form 1 — cos^ 6b for signal and cross-feed 
components, and flat for qq background. 



C. Fit Results 

Table UTTl shows the results of the fits to the CLEO III data. All errors shown are statistical 
only. The apparently large yields of qq background reflect the large background-normalizing 
sidebands in Mb and AE and are not indicative of S/B. Typically S/B 1 for the observed 
modes. 



VI. SYSTEMATIC UNCERTAINTIES 

The net uncertainty in our branching ratio determinations is dominated by the statistical 
errors in the event yields but also includes a systematic contribution. 

We categorize systematic uncertainties in two groups, multiplicative and additive. Addi- 
tive uncertainties are those that affect the overall yield of signal events, while multiplicative 
are those that enter as scale factors in converting the yield to a branching ratio. In view of 
the following equation, 

Tyobserved 

Nbb X (eff) X (secondary BR) ^ ^ 

the multiplicative uncertainties correspond to the uncertainty in our knowledge of the abso- 
lute number of BB pairs in the data sample, denoted Nbs, and the reconstruction efficiency 
of each mode. In practice the uncertainties in the secondary branching ratios of vr'' 'j'j, 
Kg 7r"*"7r~, A — > pvr and —>■ [Kn, Ktctc^, Kmnr) are negligibly small compared 
to uncertainties in Nbb a^id reconstruction efficiency. 
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TABLE II: Functional forms used in likelihood fits. See text for discussion of terms. Linear 
combinations are indicated by coefficients a and b. 



Mode 


Fit Component M{Mb) 


t (A£/j 






Signal 


G 


G 


G 




qq 


A 


L 


Fo 




Cross-feed 


G 


G 


G 


1 n 


Signal 


G 


aGi + bG2 


G 




Cross-feed 


G 


G 


G 




Signal 


G 


aGi + bG2 


G 




Cross-feed 


G 


G 


G 


K+TT- 


Signal 


G 


G 


G 




Signal 


G 


aGi + 6G2 


G 




Signal 


G 


aGi + bG2 


G 




Signal 


G 


G 


G 


K+K- 


Signal 


G 


G 


G 




Cross-feed 


G 


G 


G 


K^K~ 


Signal 


G 


aGi + 6G2 


G 




Signal 


G 


aGi + 6G2 


G 


PP 


Signal 


G 


G 


aGi + 6G2 


pA 


Signal 


G 


G 


aGi + 6G2 


AA 


Signal 


G 


G 


aGi + bG2^ 


D^TT- 


Signal 


G 


G 


G 


D^K- 


Signal 


G 


G 


G 




Cross-feed 


G 


G 


G 



"The qq PDFs are the same for all modes. For brevity we omit this line in subsequent entries. 
^Set A includes 7^. 

A. Additive Systematic Uncertainties 

The accuracy of the signal yield obtained from the likelihood fit depends primarily on 
the fidelity of the PDFs used in the fit. A secondary consideration is the correctness of the 
product form assumed in Eq. |2l which ignores any correlations among the four fit variable 
distributions. Such correlations however are expected to be small, and Monte Carlo tests of 
the fit procedure confirm this expectation. We therefore focus on the systematic uncertain- 
ties in signal yield which arise from systematic uncertainties in the PDF parametrizations 
already noted in Section lVBl To evaluate these uncertainties we refit the data multiple times 
with one PDF parameter varied each time. The resulting signal variations are summed in 
quadrature, separately for negative and positive yield variations, ignoring any correlations 
which may exist among the parameters. A representative set of these uncertainties are dis- 
played in Table IIVI for the Ktc mode; details will vary from mode to mode. (For simplicity 
of presentation we have combined results from Set A and Set B, and merged the three 
component terms of the Fisher PDF.) The essential feature, however, is that the net addi- 
tive systematic error corresponds to a relative error of 3.5% which is substantially smaller 
than that statistical error, and also smaller than the multiplicative systematic errors to be 
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TABLE III: Results of likelihood fits: raw event yields with statistical errors. A dash in the last 
column means no cross- feed term was used in the fit. 



Modp 


Set Eff (%) 


r\ 1 cm fi 1 


qq Uivg 






A 


39.0 


IT 0+5.5 
7-814.5 


1750±42 


n+4.6 
3.9_3 J 




B 


45.3 


A Q+4.1 


1955±44 


Q+3^4 
^•°-2.3 


+ 


A 


34.9 


o o~ho.o 

2-8li,9 


1158±34 


9.3 ± 7.0 




B 


37.5 


5.7±5.9 


1139±34 


0.0 lb 3.5 



vr vr 


A 


22.1 


n r,+2.5 

2-2Ii.5 


134ibl2 


3.6t|f 




B 


22.4 


0-413.4 


211ibl5 


n c+1-7 
0-5_2.6 


K TT 


A 


37.9 


no 1+6.8 








B 


45.3 


1 n 1 +5.3 
19-1-4.6 


1 /I 1 /I 

io48±43 




T^O + 

K ^TT^ 


A 


12.3 


1 o 1 +4.4 
lz.l_3 J 


398±20 






B 


12.8 


9 Q+1-8 

. y 2 Y 


395±z0 







A 


32.6 


^ „ ^+6.2 


735±27 






B 


35.3 


IO.8I4 1 


Ton 1 00 

78U±z8 




XV vr 


A 


9.6 


q c:+2.8 

"J-'J-i.g 


iD4±io 






B 


10.5 


9 Q+2-4 
^-^-1.6 


1 Q9_l_ 1 9 

iozitiz 




K 


A 


35.2 


n q+3.7 
/.5_2.9 


945±31 


rj 0+4.2 

7-213.4 




B 


42.1 


n n_i_n v 




9 n+2-6 
2-0_i.6 


K K 


A 


13.8 


0.0±i.5 


37i±19 






B 


13.0 


0.0±0.6 


369±19 






A 


8.1 


n n-i-D 1^ 

U.UltU.O 


o4lItD 






B 


8.0 


0.0±0.5 


37±6 




pp 


A 


31.5 


0.0±0.7 


38±6 






B 


34.3 




18±5 




pA 


A 


21.9 




46±7 






B 


21.6 


0.0±0.8 


44±7 




AA 


A 


14.3 


0.0±0.7 


25±5 






B 


13.4 


0.0±0.6 


29±6 





discussed next. This pattern holds true for all modes. 

B. Multiplicative Systematic Uncertainties 

We summarize the multiplicative systematics in Table IVl 

The absolute number of BB pairs in the data sample sets the scale for all branching 
ratios. We determine this number by three different methods: counting decays of the type 
B — * D^T[~ , fitting distributions of the Fox-Wolframfl^ event shape variable i?2, and direct 
computation from the run-by-run integrated luminosities, beam energies, and the shape of 
the T(4S) resonance (normalizied to 1.07 nb at the peak). The R2 method was used in 
previous CLEO II publications 4], and in principle has excellent statistical power and small 
systematic uncertainties, but requires substantial off-resonance data that was not available 
in the first 30% of the CLEO III running period. Where off-resonance data is available, the 
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TABLE IV: Additive systematic errors due to PDF variations for B K'^Tr~. Entries show 
change in efficiency-corrected signal yield (events) resulting from a parameter variation of one 
standard deviation up (high) or down (low). L and R refer to left and right sides of an asymmetric 
Gaussian distribution. 





Parameter 


Result of Parameter Variation 








T "J- * 

Low-variation 


High-variation 




Mb 


mean 


—0.1 


—0.1 






width 


— 1.4 


+1.3 


a 


A IP 


mean 


— \.6 


+1.2 


bJD 




width 


— Z.O 






XT 


mean 


-1.0 


+1.0 






width (L) 


-0.3 


+0.3 






width (R) 


-0.8 


+0.8 




Mb 


A 


-1.1 


+1.1 


a 


AE 


slope 


-0.1 


+0.1 


o 

bjO 


XT 


mean 


-0.6 


+0.6 


u 
03 




width (L) 


-0.6 


+0.6 






width (R) 


-1.0 


+1.0 






areas 


-0.8 


+0.8 




Total 


-4.1 


+3.8 



Dtt method and the i?2 method agree very well, and since the D-k method is available for 
all data sets we use it. The direct computation technique is used only as a check of the 
other methods, and is found to be in good agreement with them. In the Dtt method, three 
secondary modes are used, — > 7^~7r+, 15° K~'K~^7r'^, and K~'K~^'K~7r~^ , and a 

small cross-feed from B — > DK is subtracted. 

To avoid the additional uncertainties implied by secondary branching ratios, we em- 
ploy CLEO II Nbb determinations to set the absolute scale for CLEO III: 



n{BB)jjj ^ n{Dn)iii en 
Tl{BB)j, n(D7r)ii em 



(5) 



Event rates per unit luminosity (7^) and efficiencies (e) are determined separately for CLEO 
II (subscript //) and CLEO III (subscript ///) , and for each of the three secondary decay 
modes. In the end the dominant limiting uncertainty in this technique is the statistical error 
in Dtt yields. 

Rare B decay modes involving tt^, Kg, or A in the final state have additional uncertainties 
associated with the efficiency to reconstruct these particles. We determine the reconstruc- 
tion efficiencies in Monte Carlo (eMc) simulation and then perform a separate determination 

in data (cdata)- The total error in the ratio eDATA/cMCj which includes both statistical 
errors and some systematic errors (such as branching ratios) is then interpreted as the sys- 
tematic uncertainty in the reconstruction efficiency. For the data determination consists 
of measuring the ratio 

_ iV(L>° ^ X-7r+7r°)/g(ii--7r+7r°) 
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TABLE V: Multiplicative systematic errors. Entries show the fractional change in branching ratios 
for each contributing source. Entries above the line affect all modes while those below only affect 
modes involving the corresponding particles. All values quoted are for CLEO III. 



Source of uncertainty 


AB/B (%) 


Absolute number of BB pairs 


8% 


Monte Carlo Statistics 


1% 


Single Track Reconstruction Efficiency 


1% 


Particle ID Efficiency per Identified Track 


3% 


Single vr'^ Reconstruction Efficiency 


10% 


Single Kg Reconstruction Efficiency 


7% 


Single A Reconstruction Efficiency 


17% 



where we take the ratio of branching ratios obtained from reference [9|] to be 3.44 ± 
0.22. We find eoATA/eMC = 1-00 ± 0.08 ± 0.02 where the second error reflects conservative 
uncertainty in the Dalitz amplitudes of —>■ K'tt'^tt^. A similar study was done using 
7] — *• 77, rj —>■ 7r°7r°7r°, and 7] —>■ 7T^7r~n^ decays. We anticipate that further study will refine 
the 7r° systematic error estimates. A more precise determination of the systematic error, 
however, is not called for by this analysis as any uncertainty under ~ 20% changes our 
signal sensitivities only marginally. Kg reconstruction uncertainty is determined similarly 
from comparing — >• K^n^ and — > K'n'^n'^, which yields eDATA/cMC = 1-01 ± 0.07. 
In the case of A the comparison is of A^ — > Avr and Ac pKn, and we obtain eDATA/cMC = 
0.93 ± 0.17. In the A case, the net uncertainty is dominated by the relatively poorly known 
branching ratios. In all cases the systematic uncertainties estimated by this technique are 
conservative (large) but still do not dominate the final total error. 



VII. CLEO III RESULTS 

Event yields for the CLEO III data subsets A and B are given above in table IIIll Because 
the signal efficiencies of Set A and Set B differ slightly the event yields in the two datasets do 
not have exactly the same meaning and are not directly comparable or summable. To obtain 
overall CLEO III results we express the measurements of Set A and Set B in the common 
language of branching ratios, B = Hgig/ {N^j^e) forming the joint likelihood £stat('B) = 
Ca{B)Cb{B). The subscript "stat" emphasizes that this version of the likelihood function 
reflects only statistical features of the data. We fold in systematic errors, which are common 
to Set A and Set B, by convolving the normalized statistical likelihood function 

^ _ -^stat(^) 

with both additive event yield uncertainties r], distributed according to an asymmetric Gaus- 
sian, ^(77), and multiplicative scale factor uncertainties p, distributed according to a symmet- 
ric Gaussian G{p). The widths of these distributions have been discussed above. Formally 
this convolution may be written 
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For convenience the double convolution is performed by a Monte Carlo method. 

The resulting distribution of C{B) is the final CLEO III likelihood function including 
all of the uncertainties in the measurement. From it we find the minimum of the — 21n£ 
distribution to measure our mean, and find the la intersections to determine the errors. 
Since this is the total error, we unfold the systematic error by subtracting the statistical 
error in quadrature from the total error. We set 90% confidence level upper limits by 
determining the value of B for which 

rB^ 

/ C(p)dp = 0.90, 
Jo 

and calculate significances by looking at the zero yield value of the — 21n£ distribution. In 
the limit of a purely Gaussian likelihood function, this definition of significance reduces to 
the signal yield divided by its one standard deviation error. 



VIII. COMBINED CLEO II AND CLEO III RESULTS 

We combine CLEO II and CLEO III measurements using the likelihood functions de- 
scribed above for CLEO III and reported in Ref. jT^ for CLEO II. For some modes we use 
previously unpublished likelihood functions. The baryonic modes and the K^K'^ mode were 
analyzed here with the full CLEO II data set for the first time. 

Particle identification in CLEO II was limited, and modes with potential for K/tt misiden- 
tification, such as B ^ tt+tt^ and B K^n^ , were analyzed in terms of two-dimensional 
likelihood functions, C{Nt^+-„-, Nx+n-)- The improved K/n separation in CLEO III however 
permits us to treat these modes independently in the new data. To combine CLEO II and 
CLEO III likelihood functions, therefore, we first project the two-dimensional CLEO II func- 
tions on to one-dimensional versions, using C{x) = J C{x,y)dy, and then express in terms 
of branching ratios, C{B) = C^Nsig/i^N^ge)). Systematic errors are included following the 
same method as described for CLEO III above to obtain a total CLEO II likelihood function 
for each mode. A final combined CLEO II and CLEO III likelihood function is then formed 
from the joint likelihood, /^finail-B) = 'CcLEon{B)CchEoui{B). For vrvr and Kn modes with 
non-zero yields we plot the negative log-likelihood functions in Fig. ^ Likelihood functions 
for the di-baryonic modes are shown in Fig. El Table IVII summarizes the final results, with 
separate entries for CLEO II results (extracted from the references and reproduced here for 
the convenience of the reader), CLEO III results, and the combined CLEO II and CLEO III 
results. 



IX. PHYSICALLY INTERESTING RATIOS AND THE PHASE OF Kb 

As discussed in the introduction, it is possible to extract information about the phase of 
Vub from these charmless B decay data. The method of Reference 3 is based on two ratios 
of the branching fractions which we have measured and reported above. Using the notation 
of this reference, and combining statistical and systematic errors in quadrature, the ratios 
are found to be: 

B(B^ K^Tc^) 
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FIG. 1: -21n(£/£ 

max) distributions for CLEO II and CLEO III combined, for Ktt and vrvr modes 
with non-zero yield. 

and 

1/2 

= 0.18 ±0.04. (9) 

We see that the precision available with the CLEO data is about 20-30% in these quan- 
tities. With data from the BABAR and Belle experiment s(l| we can make world (weighted) 
averages of branching ratios and reach 10-15% experimental precision in the critical ratios: 
-R* = 0.71 ± 0.09 and eexp = 0.21 ± 0.02. These numbers in turn indicate a preferred re- 
gion for 7 = Arg{Vub) which is greater than 90° Using these world-averaged data we 
construct contours in the p — r] plane according to prescription of Reference 3 and display 
the result in Figure El The dark band represents the experimental central value convolved 
with theoretical uncertainties; lighter bands show the additional coverage when 68% and 
95% experimental confidence regions are included. For reference we also overlay 68% and 
95% confidence level ellipses of the preferred apex of the unitarity triangle as obtained in 
a standard analysis based on B mixing, sin 2/3, Vub, and kaon decays|17|]. An intriguing 
discrepancy between these regions is noticeable. In the short term the most substantial 
progress to be made will be in reducing the statistical errors on the branching ratios of 
charmless B decay modes. If discrepancies survive there could be non-trivial implications 
for the Standard Model, as discussed in Ref. 3. 



^exp = tan 9c 



fx 

u 



X. THE B{B DK)/B{B Dtt) RATIO 

In view of the good K/tt separation in CLEO III data we also report a new determina- 
tion of the ratio B{B^ D^K^)/B{B^ D^n^) which benefits substantially from good 
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FIG. 2: Likelihood functions for B pp, B pA, and B AA. 



TABLE VL Experimental results for CLEO II, CLEO III, and both datasets combined. Signifi- 
cances include systematic errors. Note that the pp analysis in Ref. 4 was done in only a subset 
of the full CLEO II dataset, so the "combined" result is simply the CLEO III upper limit. Upper 
limits are 90% confidence level. CLEO II results are taken from Ref. 4, except for the K^K^, pA 
and AA modes which were analyzed in this work with the full CLEO II dataset for the first time. 





CLEO II - Ref. 4 


CLEO III 


Combined 


Mode 


Significance ^3x10^ 


Significance BxlO^ 


Significance BxlO^ 


vr+vrO 


4 9 A Q+l.6+0.5 
"i-J-l. 4-0.5 
09 c (^+2.6+1.7 
0.^ <J.U_2 3_i 7 

2.0 (< 5.7) 


2.6 4.81-1- 
9 1 H+2.8+0.8 

Z.i 'J-^-2.0-0.3 

1.8 (< 7.6) 


A A A r;+1.4+U.5 

4.4 4. 0_]^ 2-0.4 

2.5 (< 4.4) 


K+-K- 


19 17 9+2.5+1.2 

iZ ±(.z_2 4_]^2 

7.6 18.211^+1.6 

6.1 11.6+^-?+li 
4.9 14.6l|?l|^ 


> 7 19.51-1- 
4.6 20.5+^-^+1? 
5.0 13.5lt°l?:^ 
3.8 11.0l|6±2.5 


> 7 18.01-ii:^ 

> 7 18.811^12:1 

> 7 12.9+lt+li 

5.0 i2.8ll°ll:I 


K+K- 


(< 1-9) 
(<5.1) 

(<6.1) 


(< 3.0) 
(< 5.0) 
(< 5.2) 


(< 0.8) 
(< 3.3) 
(< 3.3) 


pp 

pA 
AA 


(< 7.0) 
(< 2.0) 
(< 1-8) 


(< 1-4) 
(< 3.2) 
(< 4.2) 


« 1.4) 
(< 1.5) 
(< 1.2) 
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FIG. 3: Confidence contours in the rj — p plane. The shaded bands represent regions allowed 
by the world-averaged charmless B decay measurements while the ellipsoids represent 68% and 
95% contours from conventional global fits to heavy quark measurements. The dark shaded re- 
gion corresponds to the experimental central values of the charmless data, smeared by theoretical 
uncertainty. See text for details and references. 



particle identification. The original CLEO II publication is available in Ref. pi 

For this analysis, candidates are reconstructed in three secondary modes, — >■ 
i^~7r+, — >■ K'TT^n^, and — > K~n^7c~7T^ . Requirements for the B — > D^h^ modes 
include a 30 MeV mass cut, a 100 MeV AE cut, and standard particle ID as described 
previously on both the primary h~ from the B and on the secondary kaon from the D^. The 



7r° mass for the D^—^K' 



mode is required to be within 30 MeV of its nominal mass. 



For the D^K likelihood fit, D^n is included as a cross-feed background, and corresponds to 
approximately 50% of the DK yield shown in Fig. |3] (Both D*7t and Dp were found not 
to be significant backgrounds to either signal mode.) Fig. 0] shows Mb distributions for 
B^ D^7r~ and B~ — > D^K^ candidates with the likelihood fit shape superimposed. 
Combining the three submodes, we find 



B{B-^D^K' 
B{B-^D^7i- 



(9.9: 



-1.4- 
-1.2- 



-0.7\ 
0.6) 



XlO' 



(10) 



Most systematic errors cancel in this ratio, with only a small residual arising from the 
particle identification requirements imposed on the primary n/K in both numerator and 
denominator. 



XI. SUMMARY 

We have presented final results from the CLEO experiment on charmless hadronic B 
decays. The decay modes include the ten rni, Kn, and KK final states as well as the 
dibaryonic states pp, pA and AA. In addition we have presented a new determination of the 
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FIG. 4: CLEO III data: the Mb distribution for D^n- (left) and D^R- (right) 

candidates. 

ratio of branching ratios B{B —>■ DK)/B{B Dn). The results are based on the full CLEO 
II and CLEO III data samples totalling 15.3 fb^^ at the T(AS), and supercede previously 
published results by this collaboration. 
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